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ABSTRACT 


Finite element, elastic and plastic analyses of various 


aluminum panels, containing holes and notches, were conducted 
for comparison with photoelastic experimental results. A 
FORTRAN IV program, ADINA (Automatic Dynamic Incremental 
Nonlinear Analysis), was used for both linear and nonlinear 
analyses. Mesh refinements were used for each panel and the 
monotonically convergent results were extrapolated using 
Richardson's method. Stresses were locally smoothed from 
the Gauss integration points to the nodal points. Eight 
noded, isoparametric elements were used throughout. Modi- 
fication to an ADINA preprocessor program, also coded in 
FORTRAN IV, was made for use with a VERSATEC plotter. 

Comparisons were made to the elastic, analytic series 
solution by Howland for a circular hole in a finite strip. 
The finite element results varied by less than one percent 
from Howland's solution. Handbook values for the elastic 
stress concentration factors of the geometries investigated 
differ from finite element results by less than one percent 
in all cases. The photoelastic works of Frocht were also 
compared where applicable. Stresses in the plastic range 
obtained from slip-line theory for a rigid-perfectly-plastic 
material show excellent correlation to a finite element 
analysis of such a material. Comparisons to elastic and 
plastic experimental data were made for the panels 


analyzed and show good correlation to finite element results. 
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I. INTRODUCTION 


Development of on-board fatigue monitoring systems 
for Naval aircraft have made it possible to record exten- 
sive structural loading data in flight. The strain gages 
used in such a system must be located away from stress 
concentration areas to prevent their fatigue; however, 
these areas are of the greatest interest in analyzing 
and predicting fatigue life of the structure. Understand- 
ing the relationship between nominal, far-field stresses 
and local stresses in critical areas thus becomes vitally 
important. Recent experimental investigations into the 
effect of uniform, far-field loads on stress concentration 
areas have been made at the Naval Postgraduate School 
(NPS) using photoelastic techniques [Refs. 1, 2 and 3]. 
These experiments involved loading 7075-T6 aluminum into 
both the elastic and plastic regions, as well as measure- 
ments of residual stresses resulting from plastic yielding. 

Finite element analyses (FEA) of the aluminum panels 
used in the experiments of Stenstrom [Ref. 1] were con- 
ducted. The panels used in the experiments of Engle 
[Ref. 2] and Stuart [Ref. 3] have similar geometry. The 
finite element programs available at NPS were surveyed 
and ADINA [Ref. 4] was chosen for its proven ability to 


produce the nonlinear analyses required for plastic 
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yielding of aluminum. To provide increased accuracy, 
each panel was modeled using two meshes. The results 
obtained for the coarse and fine meshes were extrapolated 
to a final result using the Richardson extrapolation 
technique [Refs. 5 and 6]. 

Along with ADINA, a preprocessor program, PSAP1 [Ref. 7], 
was used to verify mesh connectivity prior to analysis by 
ADINA. PSAP1 provides a graphical output of the finite 
element mesh and was coded for the CALCOMP plotter installed 
at NPS prior to 1978. For this thesis PSAPl was adapted 
for use with the VERSATEC plotter now installed at NPS. 

The stress-strain material properties of the 7075-T6 
aluminum actually used to make the panels had to be estab- 
lished to provide an accurate material model for use with 
ADINA. Material testing was conducted to establish the 
Young's Modulus (E), Poisson's Ratio (v), yield stress (a), 


strain hardening modulus (E,) and the Ramberg-Osgood co- 


foe anit Kacbet sence et ties poten iy. th pharm tava Rallitierenents aliens ngmatemanin ot 


efficients 8 and n. 

Comparisons were made to other works, in addition to 
the experiments conducted at NPS. The initial analysis 
involved a comparison of FEA to the results of Howland 
[Ref. 8], for a circular hole in a finite strip, to vali- 
date the methods used. A comparison of FEA to plane 
stress, Slip-line theory, for rigid-perfectly-plastic 
Material was also included as a validation for the plastic 


analyses. 
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II. MATERIAL PROPERTY TESTING OF 7075-T6 ALUMINUM 


The elastic and plastic material properties of the ’ 
aluminum panels were established by tensile tests of uni- 
axial specimens made from the same mill run. The speci- 
mens were manufactured and tested according to current 
ASTM standards [Ref. 9]. MICRO-MEASUREMENTS, EA~-13-125AD- 
120, precision strain gages with a temperature compensated 
bridge circuit were used on all specimens. Transverse 
gage sensitivity errors were corrected according to the 
manufacturer's recommendations [Ref. 10]. Critical cross- 


section measurements were made with a micrometer. 


A. TESTS FOR AXTAL LOADING 

Initial tests of the two-gaged specimen, Fig. 1, in 
the MTS testing machine indicated a significant bending 
moment was being produced by the 30,000 1b GRIFF grips. 
To investigate this problem further, tests were conducted 
on both the MTS and RIEHLE test machines with a five- 
gaged specimen shown in Fig. 2. The results of these 
axial loading tests, shown in Table I, verified that the 
GRIFF grips on the MTS test machine do not give axial 
loading. An inspectionof the gripped region on the speci- 
men indicated that the jaws of the grip were not applying 


a uniformly distributed force and thereby induced a bending 
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moment by off-axis loading as shown in Fig. 2. The grips 
on the RIEHLE test machine gripped evenly and a uniform 


strain distribution resulted as seen in Table I. 


B. CHARACTERISTICS OF 7075-T6 ALUMINUM PANELS 
The following characteristic properties were determined 
from the four specimens tested. 
1. Young's Modulus (E) 
Tests were conducted using the specimen shown in 
Fig. 3 on the MTS test machine with 10,000 1b INSTRON grips, 


which gripped the specimen evenly. The results of testing 


mee 


three specimens are shown in Tables II to IV. Linear 
regression in the elastic range of all the test results 


6 psi, with a 


determined a Young's Modulus of 10.12 x 10 
correlation coefficient of 0.9996. 
2. Poisson's Ratio (v) 


Tests were conducted using the specimen shown in 


ye Fig. 1 on the RIEHLE test machine with 10,000 1b RIEHLE 


Re ee ee a ag a ae ee Oe ee een oor, 


. grips. The results are tabulated in Table V. Linear 


regression of these results in the elastic region deter- 


mined Poisson's Ratio to be 0.3256 with a correlation if 


coefficient of 0.99996. 


3. Yield Stress and Strain Hardening Modulus 


These values, required for ADINA's bi-linear 
material model, were determined graphically using the 


data from Tables III and IV. Plastic region data in 
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Table II is not reliable because of excessive creep en- 


countered during that test. 


0.2% offset yield stress, o,, 76,000 psi 


strain hardening modulus, E. = 566,000 psi 
The graphical fit of these values to the test data can be 
seen in Fig. 4. 
4. Ramberg-Osgood Coefficients 
The Ramberg-Osgood equation for elastic-plastic 


een noni tial amnedeaibemmeicediaie, 


stress-strain characterization is given by: : 


_ oO o,n 
e=5 t+ 8(5) (1) 


where: 
€ = strain 
ag = stress 
E = Young's modulus. 
The 8 and n coefficients were determined graphically from 
) the data of Table IV, by the method given by Rivello [Ref. 11]. 
- The data in Table IV gave the following values which are 
the best fit to the combined test data 


43 


8 = 1.479 x lo fe 


a n = 21.58 


The graphical fit of these values, in Eq. (1), with the 


test data is also shown in Fig. 4. 


| 
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MODIFICATION TO GRAPHICAL PREPROCESSOR 


The use of a graphical preprocessor program, such as 


PSAP1, is vital in detecting mesh errors that may other- 


wise go unnoticed. Establishing the proper node locations 


and element geometry prior to analysis for a complex code 


such as ADINA is of utmost importance. 


A. PSAP1 MODIFICATIONS 


The program PSAP1 was originally coded in FORTRAN IV 


for use on the NPS IBM 360/370 installation with the 


CALCOMP Model 765 drum plotter. The CALCOMP system o: 


installed at NPS used the +Y axis as the unlimited plotting 


axis, see Fig. 5. The entire plotting logic in PSAP1 uses 


this orientation of axes to allow multiple plots in a con- 


tinuous strip. With the VERSATEC Model 8222A electrostatic 


plotter now installed at NPS, the +X axis becomes the un- 


limited plotting axis, shown in Fig. 5. To avoid an ex- 


tensive recoding of PSAP1 for use with the VERSATEC plotter, 


a simple coordinate transformation of the plot was made in 


a limited number of short subroutines. First, all installa- 


tion dependent plotting calls used in PSAP1 were identified. 


These involved seven plotter functions for which new sub- 


routines were coded. 


4 
: 


Function New Subroutine 


50 5 a ig lin i aa a a cada ae Se 


Initialize Plotter CALCMP 
Move Plotter Pen CALPLT 
Letter on Plot NOTATE 
Number on Plot CALNUM 
Determine Current Pen Location CALWH 
Draw a Line CALINE 
Stop Plotter PSTOP 
The subroutines listed above merely rotates the plot 
to coincide with the VERSATEC axis orientation and retain 


all features originally in PSAP1. Since all plotter 
hardware code is now isolated in these seven subroutines, 
future adaptations to other plotting systems is simplified. 
To provide documentation of this update to PSAPl, a complete 


listing of the new program is provided in Appendix D. 


B. USE OF PSAPI 

Previous use of PSAPl1 on the IBM 360 system necessitated 
use of a load module since PSAP1 took over one minute of 
CPU time to compile. With the new IBM 3033 system 
compilation requires eight seconds; however, use of a load 
module or disk stored source code is still recommended 
since PSAP1l contains roughly 2,500 lines of code. Appendix 
A contains sample JCL to use PSAP1 on the IBM 3033 MVS 


system. With minimal effort PSAP1 could also be set up 
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for use on the IBM VM/370 system. The user's manual for 


PSAP1l is in Ref. 7. In addition to a mesh plot, PSAP1 


provides a listing of the node coordinates, element 


connectivity and several key input values used in execu- 
tion of ADINA. This information provides a useful check 


of the input data. 


perme suamsneneiamneneenamenesen es . = — 


IV. FINITE ELEMENT ANALYSIS (FEA) 


| A. DESCRIPTION OF MESHES USED 


1. Length to Width Ratios and Boundary Conditions 


Initial finite element models of specimens had 


length to width ratios near one, as used by Garske [Ref. 12] 
for his FEA, but they did not provide the desired uniform 
distribution at the loading boundary. Specimen length to 
width ratios of 3-5 were used by Armen, Pifko and Levine 
{Ref. 13] in their FEA and by Stenstrom [Ref. 1] in his 
photoelastic experiments. The criteria established to 
determine uniform boundary stress distribution was uni- 
formity in nodal displacements along the loaded edge as 
discussed by Segerlind [Ref. 14]. In the models used for 
FEA in this thesis, nodal displacements were uniform to 

‘ within 0.1%, and the resulting stress distribution was 
uniform axially to within 0.1% at the panel ends. In all 
cases two-dimensional, eight noded, isoparametric elements 


were used. These higher order elements cannot be loaded 


in an "intuitive" manner as discussed by Zienkiewicz 
{Ref. 15, p. 223]. Figure 6 shows the nodal loading 


required to obtain a uniform surface load. 
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2. Element Mesnes 
Two meshes were developed for each panel analyzed. 
A reasonable effort was made to keep element corner angles 
as close to 90° as possible to reduce the effect of element 
distortion discussed by Hopkins and Gifford [Ref. 16]. 


All meshes modeled a quarter of the actual panel by using 


TS deed oe 


the two axes of symmetry as is common practice in FEA. 

The step from course to fine element meshes was made so 
that each element in the course mesh was subdivided into 
four smaller elements of the same type. Such a mesh sub- 
division can be expected to give monotonic convergence of 
results, Cook [Ref. 17], and allow extrapolation to results 
of an infinitely fine mesh. Figures 7 through 13 illus- 
trate the element meshes used in this analysis as plotted 


by PSAP1. 


B. COMPUTATIONAL PROCEDURES 

1. Using ADINA 

Once the mesh has been developed, input data is 

prepared in accordance with the ADINA user's manual [Ref. 4]. 
This same set of data is then used as input for PSAP1 to ~ 4 
check for errors and provide a graphical display of the 
element mesh. After preprocessing by PSAPI1, the data is 
entered into ADINA for analysis. In the case of linear 
analysis, two types of stress output may be specified, 


nodal point or Gauss integration point. Nodal point output 
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can be computed for up to eight node point stresses for 

each element. Since 2x2 Gauss integration was used, four 
Gauss point stresses were computed for each element. The 

2x2 Gauss integration is recognized as the most efficient 
integration order for this type of analysis [Ref. 15, p. 284]. 
The linear analysis used an isotropic linear elastic material 
model (MODEL "1" in Ref. 4) which required input of E and 

v material properties. The nonlinear analysis allows only 
Gauss point stress outputs and uses a bilinear elastic- 
plastic material model, with von Mises yield condition and 
isotropic strain hardening (MODEL "8" in Ref. 4). 

For static analyses ADINA uses a time function 
method to apply loads in steps. Linear analysis loading 
was accomplished in a single step to a nominal value of 
3,000 lbs load. Nonlinear analysis loads were applied in 
ten steps to a maximum value, matching the experimental 
loads, and then unloaded to zero in ten steps to obtain 
residual stresses. The stress output from ADINA is a 
listing of nodal or Gauss point stresses for each element. 
Since the only area of interest in this analysis was the 
distribution of stresses along the reduced cross-section, 
no large post-processing program was developed or used. 
All final computations using ADINA output data were accom- 
plished on a HEWLETT-PACKARD 9830A calculator, using short 


programs coded in BASIC. If more extensive stress 
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distribution information were desired, some form of auto- 
mated post~processing would be necessary to reduce the 
computational workload. At a minimum, nodal stress outputs 
by ADINA must be averaged to obtain unique values of stress 
at nodes shared by more than one element. 
2. Richardson Extrapolation 

The use of course and fine meshes allows extrapola- 

tion to an infinitely fine mesh as discussed earlier. 


Richardson extrapolation [Ref. 5] was used in this analysis 


where: 
a (hi)™ = oi(h_)™ 4 
g oF Zs (2) i 
extrap hm}, ™ ; 
B ¢ , 
iq 
where ‘ 
if 
= 1 iq 
) Cextrap extrapolated solution i. 
oa = solution obtained with h=h, : 
on = solution obtained with h=h, ft 
i 
: he = linear dimension of course element : 
ha = linear dimension of fine element d 
\ 
m = 2 (for this analysis) * 
The exponent m is determined by the order of the discre- : 
tization error O(h"). Since h represents the length of an : 


element the element area is represented h?, In a two 


dimensional problem such as this O(h™) is of the order of 


{ h?, the area of an element. In the mesh refinement scheme 
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h 
used h, = sh or = = 7 Equation (2) can be rewritten 
7 c 
of (SE)? =o (Pe)? 
c ne F he 
a = (3) 
extrap hp iy 
ae - (-=) 
c c 
thus lL 
Oo, - FO 
ee 3 4 “c 
Pextrap es E} 3) 
4 


Equation (4) then becomes the relation to obtain extrapolated 
stresses from coarse and fine mesh results in a two dimen- 
sional analysis. Better extrapolations can be obtained by 
using three or more refined meshes, but, the computational 
effort increases significantly. 

3. Optimal Stress Locations and Local Smoothing | 

It is generally accepted that the most accurate : 

sampling points for stresses are the Gauss integration i 
points within the element [Ref. 15, p. 281, and Ref 18]. 
In this analysis, the nodal points are of the greatest i 
interest; thus a technique of local smoothing must be ; 
applied to the integration point stresses to obtain nodal ~ oH 
stresses as reported by Hinton and Campbell [Ref. 19]. 
The formulation of this local smoothing technique for ADINA 
elements is developed in Appendix B. The nodal values 4 
obtained must then be averaged if shared by two or more 


{ elements. 
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4. Computational Times 


Because of the extremely large size of ADINA (about 
17,000 lines of code in the NPS version) and the out of core 
solver, it does not adapt well to time sharing systems. 
Using the IBM 360 system at NPS, ADINA required 31 user 
defined overlays to create a manageable load module in 
about 30 minutes of CPU time. With the new IBM 3033 MVS 
system at NPS, ADINA is compiled without overlays in about 
one minute. When using a load module, the program execu- 
tion took less than 2 minutes CPU time on the IBM 3033. 

In addition to ADINA, the preprocessor (PSAP1) 
and post-processing techniques involve considerable time 
and effort. Figure 14 is a flow chart of the computational 
procedure used in this analysis. An example of the JCL 
to use ADINA in load module form on the IBM 3033 MVS 


system with use of the mass storage facilities is shown 


in Appendix C. 
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V. RESULTS OF ANALYSIS 


A. CIRCULAR HOLES IN LINEAR MATERIAL 

The FEA results for a circular hole in a finite width 
strip were used to validate the elastic computational pro- 
cedure discussed earlier. The results of Howland [Ref. 8] 
were compared to both the Gauss point smoothed results and 
the nodal output results in Fig. 15. The stress concentra- 
tion factor o/s, is referenced to the far-field stress. 
The smoothed results give the best match to the results 
of Howland at the edge of the hole, and the only signifi- 
cant variation between the two FEA methods occurs within 
the first 0.25 inches from the edge. In order to obtain 
the 0.25 inch stress value for the coarse mesh, in the 
Gauss point smoothed result, a midside node value had to 
be obtained by the averaging method discussed in Appendix 
B. The linear distribution of smoothed stresses alcng the 
sides of the element, [Ref. 19], appears to produce a 
less accurate result in this area of extreme stress 
gradient, when compared to ADINA's nodal output result. 
This tendency was noted in all cases; however, the peak 
stress values from smoothed results consistently gave 
better correlation with other investigators [Ref. 20]. 

A circular hole with A=0.25 was also analyzed and 


compared to the experimental results of Stenstrom [Ref. 1] 


27 


along with an interpolation of Howland's results. The me 
experimental data correlates well with the FEA results; 
however, the a. experimental data shows significant varia- 
tion between 0.125 and 0.375 inches from the edge of the 


hole, as seen in Fig. 16. 


B. OPPOSITE U NOTCHES IN LINEAR MATERIAL 
The results for linear analysis are presented in non- 
dimensional stress concentration form; however, the 
normalizing stress changes. For U notches Kp is the 
stress concentration factor referenced to the theoretical, 
nominal stress (o.) in the reduced cross-section where 
os Load/Area of Reduced Cross-Section. 
1. Shallow Notch Panel 
The FEA results plotted with the experimental 
data of Stenstrom are shown in Fig. 17. Once again both 
FEA results are shown and the variation for the two methods 
occur within 0.25 inches from the notch edge; however, 
there was less variation than was seen in the circular 
hole analyses. 
For this panel the experimental data appear to be 


uniformly below the FEA results for o The ore data shows 


g 
significant variation at the 0.625 inch point but follows 
the proper trend within 0.5 inches from the notch edge. 
According to data collected by Peterson (Ref. 20], this 


notch geometry should yield a maximum Kn = 2.74. The 
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Gauss point smoothed results matched this value exactly. 


Stuart [Ref. 3] reported a K, = 2.69 for this same notch 


geometry, with a standard deviation of 0.187 in the 14 
samples he measured by photoelastic methods. Early photo- 
elastic work by Frocht [Ref. 21] determined K, = 2.7 for 
this notch geometry. The FEA results appear to be in good 


agreement with other investigators, for this notch geometry. 


2. Deep Notch Panel 
Results of the FEA for a deep notch were plotted 


with Stenstrom's experimental data in Fig. 18. The results 
of the two FEA methods again diverge within 0.5 inches from 
the notch edge. FEA stress values at the 0.25 inch point 
have spread farther apart in this case since the stress 
gradient is very severe at that point. The experimental 


data correlates well for both o, and o.) however, the maxi- 


9 


mum experimental Tg is considerably lower with a Kn = 3.83. 


Results reported by Stuart for this notch geometry was a 


Kn = 4.05 with a standard deviation of 0.219 for 14 speci- 


mens measured by photoelastic methods. Frocht [Ref. 22] 


reported a photoelastic K, = 3.9 for this notch geometry, 


T 
but concluded that the result was 5-10% low, giving a 


corrected range of K,, from 4.1 to 4.3. The Gauss smoothed 


T 


FEA result gave a K,, = 4.24 which compares well to an 


T 
empirical relation given by Peterson [Ref. 20] for 


r/d < 0.25. 


2t 


yy 2t 2t, 2 2t,3 | |. 
Kya (1 - $=) (0.78 + 2.243vt7x) 0.993 + 0.18 S - 1.06() + 1.7115) (5) 


where 
t = notch depth (3.9375) 
r = notch radius (0.625) 
d = minimum width (15.625) 
D = maximum width (23.5) 


Inserting the values above into Eq. (5) gives a K, = 4.26, 


T 
which is 1/2% above the FEA result. Other FEA results 
reported by Armen, Pifko and Levin [Ref. 13] and Griffis 
{Ref. 23], using linear strain triangle (LST) elements, 
produced Kn values within 5% of those produced by use of | 
Eq. (5). The rectangular elements used in this analysis 

are known to give better results than LST elements as 

noted by Clough [Ref. 24]. It is clear that the FEA 


results obtained for this notch are in good agreement 


with other works. 


C. OPPOSITE U NOTCHES IN NONLINEAR MATERIAL 

The analysis for loading into the plastic region of the 
7075-T6 aluminum was made using the bilinear material 
model discussed earlier. The loads used were selected to 
match those used in the experiments of Stenstrom; thus 
allowing direct comparison. The strains obtained in those 
experiments were used to solve for stresses by use of the 


Prandtl-Reuss plastic flow equations. 


eee, fag 
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1. Shallow Notch Panel 
The results of FEA for the three load cases, 60,000, 
65,000 and 70,000 lbs are presented along with the experi- 
mental results in Figs. 19 through 21. The os results 
compare well although no trend for peak 4 stress away from 
the notch edge is shown in the experimental data. In all 
cases the FEA determined the peak Tg stress to occur near 


the yield boundary, and the gradient of the o, stress to 


8 
fall off dramatically in the plastic zone. This character- 
istic behavior of the J, stress was reported by other 
investigators (Refs. 13 and 23] using FEA on 2024-T3 aluminum. 
Plane elastic-plastic stress distributions reported by 
Frocht [Ref. 25] show similar trends. The experimental 

data also shows a marked change in the gradient ofo, stress 
within the plastic region. The growth of this plastic 
region is approximated using the FEA results for this notch 
in Figs. 22 through 24. 

Experimental data for the J, stress distribution 
matches the FEA results closely except at the notch edge 
where the measured o, does not go to zero as it should. 

The characteristic peak value of d. near the plastic 
boundary as seen in the FEA results is also shown by the 
test data. 

The FEA residual stress computed upon unloading 


from the three load cases are shown in Figs. 25 through 30. 


The characteristic distributions of the Ty residual stress 


agrees with those reported by others [Refs. 25, 26 and 27]. 
The experimental residual stress distributions reported by 
Stenstrom show similar trends but significant variations 
when compared to the FEA results. 
2. Deep Notch Panel 

Three load cases were computed to plastic loading 
levels; however, only limited experimental results are 
available for this notch as seen in Fig. 31. The 30,000 lb 
load is just at the onset of yield in the notch root area. 
Limited residual experimental data [Ref. 3] was available 
for comparison in Figs. 32 and 33 which are plots of the 
residual Te and OL stress distributions as a result of the 
three loading cases. The 30,000 lb load has caused yielding 
in a small region at the root of the notch as seen in Fig. 34. 
Figures 35 and 36 illustrate that the plastic zone does not 
grow to the extent it did in the shallow notch. The stress 
gradients are very severe in the deep notch and as a result 
the 0.25 inch sampling points used in the FEA may be useful 
in only showing gross trends close to the notch edge. The 
trends appear to be much the same as in the shallow notch, 


with the peak g. and o 


3 r stresses occurring near the yield 


boundary. The yield boundaries shown in Fig. 31 are 


approximations based on qualitative analysis of the finite 


element results. The Ty experimental data for the 30,000 


lb load case correlates well with the peak stress, again 
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appearing low as it did in the linear analysis. The residual 


Oo, Stress distribution in Fig. 32 follow much the same trends 
as seen in the shallow notch but with much higher gradients 
within the first 0.5 inches from the notch. The FEA results 
for the residual Te stresses shown in Fig. 33 indicate a 
limitation in the element size used since the O. value of 
the notch edge does not return to zero as it should. Because 
of this problem the data may be questionable for showing 
proper trends in the first 0.5 inch from the notch edge. 
3. Rigid-Perfectly Plastic Panel 

A stress distribution for the theoretical material 
used in slip-line theory was desired. By using ADINA's 
bilinear material model such a material could be approxi- 
mated reasonably weli. For this analysis a Young's Modulus 


(E) of 1076 


psi was used to model perfect rigidity. The 
strain hardening modulus (E,) was set to zero to model 
perfect plasticity. Poisson's ratio (vv) was 0.4999999, 

as close to 0.5 as the computer would allow. The E = 1076 
also represents a computer limit in approaching an infinitely 
large E. Figure 37 illustrates the results obtained and 
compares them to a slip-line solution. The results are 
normalized to the arbitrary 73,000 psi yield stress used 

in the analysis. The Ty values obtained agree exactly 

with slip-line theory. Conversely the oy results do not 


reflect the same values as slip-line theory, but do show 


a similar trend. The growth of the plastic zone obtained 


is shown in Figs. 38 through 40. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The results obtained from the FEA have proven useful 
in determining the validity of experimental data gathered 
by photoelastic techniques. The FEA results varied by 
less than 1% when compared to published analytical results 
and handbook values. Experimental data from Stenstrom's 
photoelastic work correlated well with the FEA results. 
The primary exception was the residual stress experimental 
values, which varied significantly from the FEA results. 
The variation may be in transforming the residual strains 
measured photoelastically into stresses for comparison 
with the FEA results, since ADINA only provides a stress 
output. Limitations of ADINA's bilinear material model, 
initially considered severe, do not appear to have hampered 
this investigation. A possible exception is the residual 
analyses where the transition region from elastic to 
plastic strains becomes especially important. The Gauss 
point smoothed results gave the best correlations at the 
edge singularities in all cases; however, due to the 
limitations noted at 0.25 inches from the edge, the 
results at that specific point may not be as accurate for 
this method. The nodal output results gave consistently 


higher stress values at the edge singularities. Because 
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of the severe stress gradients near the deep notch analyzed 
the use of a finer element mesh near the notch would probably 
produce better results. 

The effort involved in developing two meshes such as 
those used in this thesis is considerable. An automatic 
mesh generation capability would reduce the workload and 
allow experimentation with several types of element meshes. 

ADINA proved to be a useful and powerful program, as 
expected, but something simpler and less awkward to use may 
be all that is required for two dimensional analysis. Such 
a system is already in use at NPS but does not offer non-~ 
linear capabilities. If use of ADINA is to be continued in 
this type of investigation, a post-processing program should 
be adapted. There are programs available to post-process 
ADINA data at NPS [Refs. 28 and 29] but they would require 
modifications to work with two dimensional analyses and 
the VERSATEC plotter. 

Standardized material property testing would ease the 
inevitable task of obtaining basic material properties for 
use in analysis or experiments. Some form of automatic 
data collection with use of the MTS testing machine would 


allow testing of a larger sample population and provide 


statistically more accurate information. 
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TABLE’ I. MTS AND REIHLE 5 GAGE TEST RESULTS 
All Strains are 107° in/in 


MTS Test Machine 


! LOAD STRAINS 
lbs ey e5 e3 Ey Es 


REIHLE Test Machine 


ete ete ate gate oo ee 


LOAD STRAINS 
lbs 


1564 1540 


2000 
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TABLE II. MTS SPECIMEN A TEST RESULTS 


Cross-section = 0.03975 in 


Load, lbs. Strain 10— 
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TABLE III. MTS SPECIMEN B TEST RESULTS 
: i 22 
Cross-section = 0.03975 in 
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TABLE IV. MTS SPECIMEN C TEST RESULTS 
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TABLE VI. 


TABLE VII. 


A = 0.2 HOWLANL DATA 


DISTANCE 
FROM HOLE 


X = 0.25 HOWLAND 


DISTANCE 
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in. 


TABLE VIII. A = 0.2 FEA RESULTS - NODAL 


ao, = 5000 psi 
DISTANCE 
FROM HOLE, in. 


} = 0.2 FEA RESULTS ~- GAUSS OUTPUT 


5. = 5000 psi 


DISTANCE 
FROM HOLE, 


TABLE X. d = 0.25 FEA RESULTS - NODAL OUTPUT 
DISTANCE ARBORS pst 
FROM HOLE, in. 


TABLE XI. 0.25 FEA RESULTS - GAUSS OUTPUT 


o., = 4500 psi 


DISTANCE 
FROM HOLE, in. Oo, ,psi 


TABLE XII. SHALLOW NOTCH FEA LINEAR RESULTS - NODAL 


= 4225 psi 


DISTANCE 
FROM NOTCH, in. Tp 7psi Oo. ,psi 


1.25 s134.1 iPas. 5 
1.58 $332.4 (328.7 
1S 4564.7 LSoo. 5 
2. 48 4361.1 i4e4,7 


TABLE XIII. SHALLOW NOTCH FEA LINEAR RESULTS - GAUSS 


oo 4225 psi 


DISTANCE 
FROM NOTCH, in. Jy *psSL ee ;PS1L 
os er 
4.25 Brar.d birts t 
4,58 | feaoes 1384.3 be 
a.7s Sonn, 9 ris. 
1.468 Seal. 2 trai. 4 
1.25 | fo] Pobre ee ay : 
} 
1.54 4254.6 Ledd.2 
i °s $552.8 1551.7 { 
2,40 $553.5 


TABLE XIV. DEEP NOTCH FEA LINEAR RESULTS - NODAL 


oS 4800 psi 


DISTANCE 


FROM NOTCH, in. Fy +psi Oo. ,psi 


TABLE XV. DEEP NOTCH FEA LINEAR RESULTS - GAUSS 
a 4800 psi 


i DISTANCE 


FROM NOTCH, in. eg ales Oy psi 


a. 58 Bene. a mil. 
a a.75 SaS2,2 3581.9 
ha 

1.48 tS Ee | see etwas 

isco B436.1 Sata. 7 

1.58 Si4i1.5 2642.8 

Lore $230.1 f4at.2 

2.80 4622.48 SE3E.3 
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TABLE XVI. SHALLOW NOTCH FEA NONLINEAR 60,000 lb LOAD ; 


DISTANCE NO LOAD RESIDUALS 


FROM NOTCH,in. oO. ,psi Oo. ,psi oy +psi o ;psi 


TABLE XVII. SHALLOW NOTCH FEA NONLINEAR 65,000 lb LOAD 


NO LOAD RESIDUALS 


DISTANCE 
FROM NOTCH,in. Fy spsi Oo. ,psi Ty *PpSi o. *psi 


TABLE XVIII. SHALLOW NOTCH FEA NONLINEAR 70,0001b LOAD 


NO LOAD RESIDUALS 


DISTANCE 


FROM NOTCH,in. Tp psi o. »PSi Oo, ,psi Oo. ,psi 


Lec 2 RITeS. 9 esel3.<2 Fes. d SE1E.8 
1.38 TS341.5 eLaeg it 2669. 7 38. o 
1.78 TaS5r. SHETa.S Lido, 4 Job 
2.0 Stree Le4e8 6 1445.2 2207.8 
TABLE XIX. DEEP NOTCH FEA NONLINEAR 30,000 lb LOAD 


NO LOAD RESIDUALS 


DISTANCE 
FROM NOTCH, in. Sy »psi o. *Psi ao, ,psi o.. epsi 


TABLE XX. DEEP NOTCH FEA NONLINEAR 35,000 lb LOAD 


NO LOAD RESIDUALS ' 


DISTANCE 


FROM NOTCH,in. dy +psi Oo. ,psi 0, »psi Oo. ,psi 


tee nbeent 


| 
TABLE XXI. DEEP NOTCH FEA NONLINEAR 40,000 lb LOAD 
NO LOAD RESIDUALS 
iY DISTANCE 

‘ FROM NOTCH,in. J, spsi oD psi Cy psi J. Psi 
a,c 378.5 243335 a ec Pe | -2152.3 
i, 34 SeSr 2.8 Sesis. 4 ~e Lia. 4 2377.4 

: os S384. 8 24873.1 L468.2 143.1 
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Loe Br164 5% Bost, A aoe 7 ; 
a a 7 ¢ 

1.58 S425, $25.5 ria@.a { 

sere) S25 12.8 teerz.z ced re red. 3 

2 Ne S18re3, 3 153082,3 282.8 $78.5 
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TABLE XXII. RIGID - PERFECTLY - PLASTIC RESULTS 


ey = 73,000 psi 


DISTANCE 
FROM NOTCH,in. oO, ,psi o. ,psi 


TABLE XXIII. EXPERIMENTAL DATA ) = 0.25 HOLE LINEAR LOADING 


o, = 10,749 psi 


DISTANCE 
FROM HOLE, in. Oo, »psi a: +pSi 
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; TABLE XXIV. EXPERIMENTAL DATA SHALLOW NOTCH LINEAR LOADING 
| 
15,000 lb Load 
1 
DISTANCE 
FROM HOLE, in. a i Os +psi 
A. aS e 
| @.12k a2. 8 
‘ A, 258 PIS 7 
igsiec ase 
A, See Sa, S 
W.625 {ar9. 4 
TABLE XXV. EXPERIMENTAL DATA SHALLOW NOTCH 60,000 lb LOAD 
i NO LOAD RESIDUALS 
} DISTANCE 
: FROM NOTCH,in. oy +psi a. psi C5 ;psi o. +psi 
= : 
7 
1. ee Sride.d L3e2rs.4 
{ 
j 91 
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Wire _— pn Noe | em ER? 


TABLE XXVI. EXPERIMENTAL DATA SHALLOW NOTCH 65,000 1b LOAD 


NO LOAD RESIDUALS 


DISTANCE . ; : 
a. ,psL Ty *pSl oa +pSi 


FROM NOTCH ,in. 9 «psi : 
a 


TABLE XXVII. EXPERIMENTAL DATA SHALLOW NOTCH 70,000 1b LOAD 


j NO LOAD RESIDUALS 


DISTANCE 
. FROM NOTCH ,in. Ty *psi Oo ,psi Ty psi ox “psi 


a.258 Bd255, 3 12712.2 ~ $A T 48 ~a7S 8 

1 
i a.sr5 | SZS4e.2 L2S15.6 ~ 4d a “Latte, él 
a.coo | s25e4.6 ; 16814. -3O8T La -L 18.8 
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TABLE XXVIII. EXPERIMENTAL DATA DEEP NOTCH LINEAR LOADING 


15,000 1b Load 


DISTANCE 
FROM NOTCH,in. 


TABLE XXIX. EXPERIMENTAL DATA DEEP NOTCH 30,000 lb LOAD 


Elastic - Plastic 


DISTANCE 
FROM NOTCH,in. 
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APPENDIX B 


s LOCAL LEAST SQUARES SMOOTHING 


® corner Nodes 
© Midside Nodes 
+ 2x2 Gauss Points 


Two-Dimensional Isoparametric Element from ADINA [Ref. 4] 


The local smoothing expression from Hinton and Campbell 


[Ref. 19] in ADINA coordinates becomes 


oy Cc B B A Ty 
oP) B A Cc B Ors 
~ = x 
63 A B B Cc Srrr 
oy B Cc A B ony 
_ V3 ee | = 43 
where A= 1 + zs B=- 5 and C=1- > 


With Tye Jos 33 and om representing the smoothed corner node 


rr Srrye and Try as the unsmoothed 


stresses at the Gauss integration points, this expression 


stresses and <7., 


can be written in an equivalent form. 


a A BO B c. 

7 A BOC : Peay 
ont . Cc B A B Sry 
om B Cc B A Cry 
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The midside node stress values may be obtained by averaging 
the values at the associated corner nodes, since the distri- 
bution of the smoothed stresses is linear along the sides of 
the element. Smoothed stress values obtained by this least 
squares method should subsequently be averaged to obtain 


unique values at nodal points shared by adjacent elements. 


eer ae ea ne OD 


APPENDIX C 
ADINA JCL 


+ # <M H 
rs ee —~— = — me Re KR ee me ~ UU 4e 
od Fe ee _~ ~ % 
* Set ttt tt ot ot al ax 
& & o ~ o o o Cs CS o rs rs ro oe 
t HAN AN FIN eff: GIN AN FO offs HIN CN AN an be 4h 
+ FR se a ee ee ns es nes ow am RHH Wt 
* HU et ol el aL ML! ol ott et) MIL! ofl! o& Ulie HHHE es] 
Hq HR SRF Ue IR SR SP I ke oI HHH NF 
#2 HW>W>W>W>UI>Ww>W> >> >> © W> ££ OF 
to # AVAO FDO IOJAOSDOSDOSIOIO SOIR HO ##H LH 
+O FES er ee LU LL er ee Le LU ee LS Om Ft HH ax 
tat HFOUNUOUTOUCOUVOWVAUKAIOnagnovtonvwZot HHH OF 
Hoe 3% el MLW MU MU et OLE ML elt) eU etl stitnee ott 4 4 *% 
*#O HH BUBOZBO TPVUFOUTB EV FO ZFOBZFVUZBOwW eB HHH MH 
1 Ban) FEW TU eT LU eQWU WU QU LW TUIIO oI HHH ZH 
’ fo) FAFAZFAAFAZFAZFAZAFAFAZFAAFMN UN eFPOAO HHH 
] —o me ee) er Ye ee Ye te er YD mee IY wer HH 
’ HW HUW et oN eft af at oN af et off of tut tt ~ Rut oe 
' “Ne NOR RH 0 em mH KH AMR eH ~ 
] ~Y) —NOVUIODNONONODNOMNONONOVINMNO IlWMO ¥ wet Cw 
' ~ ” —OMO™MOMOMOMO“OMOMOMOMSC JOO re TF ~ 
Av NI = DODOANIANMACHOOOOCODOOCAOMACOAQH # HH F- fon) 
ala i] on Pet OF OH Bet met Ot Ot ot oe OH OD ot RTH OO 
co oO wz CTU NTI Aen aeppgcnracgnaeniacn ge Ket a 
ao io an AOQOWOWOWOWOWOWOWOWOWOWOQWOYOQAW ##at <q 
on zr On UOMNYNONONNNMNONINNONNNMNIO TMA #£O 1S) 
t> 0 - DD et tt ee SN a 
nN7s av VMMNNUNNNNNNINMANNNNNNNNMNWM EM HN ert Ney 
ons: —-_ = Wr ~ UNM EVYUYVUYNYIMHVvEYUEVrNrnvtoanyvy eee) tu 
Ie © © BRYN F Fee SF Ub ob SF SR I IE SEO OR HER 
OO Ft a = A OOO MD COMM HMM HOM eM HOH _ 
Ou: +I OO x FO apt FH ee OP OF OF OF OF OF OD OF eoZ eOFPOZF & RFEAHKNOT 
aera ¢ Oo a HY OF = mMOINDMNDINIDNDYVIONIBDNDHDMWIMDMND gas % evar 
aq@eorun + aati oOo OM OD ODO OM MD NM OM w He H TT 
ao Le HOw wo ~ DoSoSosSasasosososOsasasoso> #CHOOW 
> jhe Ol) ~O1? DQ AN AOUOUN#O AN KOKO HOA UTA KOUNUOANXOUOWO nF # Ue a 
QOZOKrF Ow NAD =Z=rzt XT FF FF FEZ eFYe FF *HtIVY 
b— Mee tee Ox MOV NZ RA. UL UL tL Lt LU SUL et ee LOR WHS H 
YrarZodazqw HF ZMH QUOCVQVOVOVOVYVODVOVUOVU OV0V 0G QUOK H# ITH 
Carr Ov HOM TSH OWOWOWOWOWOWOWOWOWOWOWOCOW # # JI 
1OMNIIET YO, HOWE UMLYLAUL YUL UYU LO LOU ZH HR HOOH 
(Uw>ry-adSexKyYOORYI>O Hewett NS Ne Pe OS OD Ow He NEM DOR et % 
i] NAZWA Teh FHAIADHEKOUADUHOROUONOKONANAHANANOLNNOKRBH ED 
10 © emsCYCoFOMWH © CLIK FOF OF OF O- Or MEMEOEOF- OE OK UE O>-HRHH I EK 
{wer # QOUOF # LOILOLULOLOLOLOLOLULOLOL ZLOUIE HE FE 
Mare HWWAZH OD eM eA 0O ¢A eM #O 20D eM 0D 0D 0o— #D OHHH * 
wood HVYYewWHO OO CIOdO0 0.07200 AO DO 200 OXF * 
Lu * J #0 Oo O86 8 HO OO O HO Oo O OO DD OO OFF 
oN S : id FANN AS VHRWKVRAHAWWAAAARW SAA WAAAANSE HH % 
wns tN FANN AHKAVLARVARWVB®AHARWHKVWRBwBWeaqnnnwnrwv# * 


Pa. < ae = re ~ 
> , e Si : s ecemmonnrte — 
i Sat ati oe : ae a ere 


fi-donag 


—- 


ie eats 26 


APPENDIX D 


PSAP1 LISTING 


(slololelelelelelolelolalololalolelelelelelelololelelelulelolclelelelelalelelcluielolalelel ola) 
MANOA SR OM OST OMNMASINOE OF OMNMOSNOE OC ONO SNOPMOROANM PAC © 
DOOCOCOO OCG qh eet test eat at ot tt AE INI NIN AIN NAN MO MO MOM AM MOP FFP te we 
PFOOCOCOOSCOOOS COO OSOCOS COCO GoCOooooooooooooeoocooceooeoo 
Lappe pas pam pam Jaw Jaw ue ue owe fe Jw p> Joe we oe p> pee ow Sow ew > os owe us ewe no > ews ne Po bm po po ee ee eo bo ee 
LOVOQOO OOOO DOO CU VO VOOOVOOGOO OOUO OO OU VOD OUDUOLYVOU 
BKOCCcocoocooocococoeoocoacococoooececeooooocoocecead 
POSSE OSSOESESSSOSSSOSSESSOSSO™NASOHSOSOHSaaHS06608S 


zor ° moO ~ 
or = S— ~ ~ 
= < =a be ~~ - 
| ad ond fe) [= 4 e e e po] [- > « 
aNd wo Ww =z 2 2 a <q “ear 
Zad (a) Du oO Oo =) 2 = a2« 
(4 Loos <x a eet oom — oe ec - C 
w Muww za eq - bh - en a) 
” i = < >a 1S) 1S) (=) we = uw i 
= 2we ww us wi Ww uw (=) zxaccood2 
< ow aor oO w ad w m wu w 
x hm 0 ad = Come Ome Om on] qq Om OO £ 
Ww Ww 2 o> Za 20 26 jes) a . OC WwW 
° 4 = > we cr O11 oO1 = t= ZF ww 
~” =-—<r~ @ a be > A w= « 
e Loh? | w pam rs | aad - Lad WwW Se wo w wl 
= 2) bn af oa <a V2 VZ us ZT eteNZTE Fa 
a ecu tw ws UJeme Lhe UL eet mwe— wn 
« - 4 row “zo (oa) od a 4 w Oo Ooftri O& m 
(=) <q a Ou [s) we 4) tt) tas? Pd ows z= 0a 
1S) (=) fie he ox ax OF Ce we OWOOE mw 
paw “ere "OZ = $12 12 12 (=) O a Qa 2 
< be ub <—w Sw NW 2 ' (oe ZS um 
> é ~>MuM) “wx = = OQ zx 20 C¢ 
a C4 [TEx 17) tu, Ly | al oe aw aw aw Q W2ZWwM a 8 
[=] oce > z mi mm YW 2eezwemn 
(=) lt ey +t moO <q <q g UJ mOQmubrerr Ww uw 
uu - MULL Wier t fon | qi dw da CC edie Die YF 
lod > [4-0] >> § a Ka FO -O (4 >a [4 
uw a uw ww Q I" MQM an QO ocee2zrt 2 Oe 
md z K-or Ow WwW adc Cm Qn Onn O eden UO eX 
Qa one W2h eo A =D © an) a COweo WwW w 
Oo << “~> «eo -_ -_ asd Oo 3Se@s>“~0 46 mM 
= u [- a Me ot Ol zo 20 20 zn nwo 2 <> 
[=] <we We oO We WwW2 w2Z C4 Ow> «<« & 
” xm>0 20 2 ue == = = aKade qw 
< 2 z wocCuw aQO = SO wo Wo wa uw WwW WwW w Or 
j=] So a e < OS OD op Zz —Omda 2 Land 
om ot <r uw FQ W wer #24 fy ty we he eee 
—_ ~- Lad — dS = XX Wet £#wIO WH WO H- aaadu &—- Z> 
@ < a LO ee Net << or armrarrare DBD ACA #£D om DO 
on Lael one UL SE ed 2 ZK DQNeat Ne teNet CO DAD N*OSLE 
od = [74 = mae WY z Ou eal oad ae NWHSQ fce3S50 
uw Oo SZMQ a0 w LY NEPOoOoqCondq@o De > we oaatcow 
a“ = (2) Zum Zz eE D> wi h be - Seewe w DOE 
< S Ww <ewo 2=z iw Ow OHO KHOTINICSTH Herwowacin aa 
=z (©) Ga pat Ge) Oo a eae Za 2a 20 Ren Ne ->a 
oO Wu WO > LW om b= - Fe NWBHDNOHNZvNA 
(=) ot fv? SLVWIASCALCCANLLKS WE TF ZeHiwWw 7 
Ine ' O —OH-DBOMDOO TCAD Oe OertODmZaen 
~ aw 2 we SaeuuCuuidbur dtd Lwoweecw 4d€L.we 
a oorTt z wt uw Te uu OU CML UL COLL UL f x 
< CeUr- jas} Zz WOCWOsWOmwoe aw wood | 
” post (=) Wale ND OCOWeadew 
a. ! - -] wn owe we OW 
~ a a % + * * “On wn 2AM 
a ire [= Q Oo & ta #0 #Q & tot 
< Zz a uw 72) Ww 7, ) zzz gad 
” os < - ui Cad tot toms =) hw 
a oad (S) —” w & [=) a a ad a aad 
> oe ta = > =z wced aca 
°o WJ a - 2 2 2 2 Maw wu 
(4 a) ud 
Oo = =z 
] = < 
72) ad 2 


VOQOVOOVOVOUOOVO UDO VU OV OOOO QQ OUND OO OU OO OU OO OULU 


98 = 


SCOCOCOSOCOCOCOOCCOCCSCeCCSOCoCeesacooesoocoocooeoooceoece ceo 


ooococcoooecoocoocescoceeoocoood 
(alalololalalalalalal= | elalalololalalalalal al al alalelal—=) 


DATA 


ALLOW TO DEFAULT 


ALLOW DEFAULT 
TYPE OF OUTPUT DEVICE TO BE USED. 


q 
- 
< 
Qa 
~- 
2 
Ww 
= 
ts 
2 
< 
J 
a 
i? a) 
het 
(=) 
2 
— 
Q 
< 
Ww 
[. 4 
jan} 
—_ 
w 
az 
Ld 
P- 4 


ARCH CENTER USE ONLY 


aS OY te ara 
AN £2S#wWMNAEN 
° YL wd ee ee ek 
AruUDe# ZOSAO 2OWO OF 
= Wether Olsen come LL 


AT NPS 
*% 


0 ** 
GEOM2 WHICH ITS PREPARED BY THE USER TO 


a (s) 

HU Zed el Me NO 
7) cq WJZ ZYrMN 
eS NOade RPOOAWE Sk SCS ve 
a AW ws el Sb I Oe Od 
Se Server Swe EDaodsoawawsd 
LOPCLMOCCO dura udquowAaAd 

1 Pema eeWu Aan 

i we wa 2 aw WoO ot 

0 1 OQWaINnmrQan a (a) 

t 


a 
HROHWH HH ee 
+ Wt * 


1 
INE GEOML WHICH READS ADINA GEOMETRY 


9 
CALL SUBROJTINE GEOM9 WHICH READS SAP IV GEGMETRY DATA. 


G ON THE VALUE OF KG 
uT 


TRY 1S NOW INPUT _IN 


*x 
IDCASE = 


NVALUS 
USE ITF KGEOM = 


QQ UN VOU COVE OU OU OO QO VU OUOUVOOUV VOVUOUNVOOUNOU 


99 


Rae Ae 


lol alelelelololelelelelelelolelalelolelelclelaleleiolslolelelaleleleololalelololelelololelelelelo) 
PDKOCHMNUMLSWOM OS OMNMSE OE OTOaINMLPNONOH OMS NUM TNOF OS Canim. 
POS ROC SCO COCO GOS ede ets ot ad tt td MIRIAM AIO HMA MM MAM MO Soh Ft St 
DOO ee ered ord 0 et pred erat grand Od cart red cad cad ened od goed gcd 0nd goed at med god ood med ed ed guret pend od pong ed somnd gue oad ed guest rend oan end md pwnd ond ud cat ome 
pam pow pew pee hace nw hoe fue jum > ne ue jaw no ue ue fe fw om o> ow Joss Jue o> nw Jew» Jw = ee ees ow ews ow ow Jee oe em ow ws em» ow o> ew es ow ow ew) 
BO VOU VOU ODOOOU HO OW OUDL UU OUD OW OUND OULU OO OUOLUOOULDU 
cocococcooocoocoeoccoccocoocomocecococoecocococacescceoo 
ANDOAGOGOOOHGHGHOAODSG OOCOOASHOOGSAGGOHOGOEHAOGACHOODGA0GG0 


e 
2 
[=] 
feng 
e) ao 
za Wa, 
oun zo 
mo ES [-4 
ox cre 
OWik UL 
Ccza> me 
[% jus To) Daal [4 
aw ive 
2ike m= ") 
ee zx = 
Q. ad o2 
[ao eT. | pe) 
Weg az = 
m = Ou - 
“wane u 
ea” ON a > 
(Sy.'4- 4-4 ww 4) 
Wiese xm 
Awa ae od [=n] 
MUA bk ou Ww 
Oar eas) [-4 
WY Lua Qu: < 
=e a a. 
aru wn Ww 
omzOou z [4 
" 2a Od a. 
WY? ol — 
NZI Pt 4 [7,) 
fiom al an = 
Octmrm 4 
Om 2zO bh = 
— Su oo 
oOoeow an -_ 
LUO> 2 zr 
=~ $104 wo) =z 
Ono po} 
«© OMS 1 | ~t 
Oo wen oad <q 
ee FUSS 2> = 
are OD C4 
QO =—nmod “nw (=) 
= - Lo 
2 Or2z2 bno Ww 
Q be poe ee oO — 7 
— [. 4 us - 
— mw aZo. -& th 
< te O = =] 
OO O2—M2 oOo azo 
mm ertuikeO az Loh -4 
uw «nd - <a 
om (US oe an» 
= [ota we wn" 
Zz w“nmaqgw mu 
w =m LA a. bbe end 
Qo ru246 Ow at od 
pre fem met oes -OQ C4 
We 
WW < ” 
" oad 32 
4 < 
Oo [=] 


IS PREPARED BY THE USER 
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IN DEGREES 
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AXIS, 


MODEL ABOJT ITS X-AXIS, 
OF MODEL ABOUT ITS Y- 


ECOND). 


3=2. 
IGNATING VERTICAL AXIS OF VIEWING PLANE 


F MODEL ABNUT ITS Z-AXIS, 


2=V> 
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